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ABSTRACT 


This report examines alternative letdown guidance strategies for 
landing of a VTOL aircraft onboard a small aviation ship under adverse 
environmental conditions. Off-line computer simulation of shipboard 
landing task is utilized for assessing the relative merits of the 
proposed guidance schemes. A sum of seventy sinusoids representation 
is used to model the ship motion time histories. The touchdown 
performance of a nominal constant-rate-of-descent (CROD) letdown 
strategy serves as a benchmark for ranking the performance of the 
alternative letdown schemes. 

Analysis of ship motion time histories indicates the existence of 
an alternating sequence of quiescent and rough motions called "lulls" 
and "swells". A real time algorithm for lull/swell classification 
based upon ship motion pattern features is developed. The classifica- 
tion algorithm is used to command a go/no go signal to indicate the 
initiation and termination of an acceptable landing window. Simulation 
results show that such a go/no go pattern based letdown guidance 

strategy provides an order of magnitude improvement in touchdown perfor- 
mance over that obtained with the nominal CROD letdown scheme. 
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ORIGINAL PAGE S3 
1. INTRODUCTION OF POOR QUALITY 


An important problem in naval aviation is the landing of a VTOL 
aircraft onboard a small ship under adverse weather conditions. 
Typically, the shipboard landing problem is divided into two distinct 
phases: (1) a decelerating approach to a stationkeeping (hover) point 
near the stem, followed by (2) a letdown or descent from hover to a 
landing on the ship’s deck. A sketch illustrating the relative geometry 
of the shipboard landing scenario is shown in Figure 1. Both the 
aircraft and the ship landing pad motions are defined with respect to 
a local level reference frame (X^, Y^, Z^) centered at the average 
position of the bullseye with the three axes along the ship's stem, 
starboard and up directions, respectively. The stationkeeping point 
(x^> y^, z^) represents the position of the aircraft upon termination 
of the approach phase and prior to the initiation of the letdown 
maneuver. 



Figure 1. Shipboard Landing Scenario 



The primary task during letdown is to bring the aircraft from the 

hover point to a touchdown onboard the ship's landing pad (approximately 

2 

40 x 40 ft ) in a safe and efficient manner. As a minimum, the 
aircraft should be able to land on the deck without damaging its landing 
gear or skidding off the deck. This requirement translates into specific 
limits upon the acceptable relative translational and angular positions 
and velocities of the aircraft with respect to the landing pad at touch- 
down. Additional constraints on the ship landing pad position and 
orientation at touchdown ma be necessary from an operational viewpoint. 

The nature of these constraints would depend upon a number of factors, 
including: (a) the aircraft-control system capabilities/limitations, 

(b) landing gear specifications and post-touchdown aircraft recovery 
procedures e.g., existence of clasping or clamping mechanisms onboard 
the ship, (c) pi lot- centered human factors guidelines, and (d) environ- 
mental and safety considerations. 

Thus two questions are posed: 

1. What is an acceptable landing? 

2. What percentage of the touchdowns would be acceptable if 
only a constant rate-of-descent letdown guidance is 
commanded to the aircraft? 

The first question may be answered by introducing the concept of an 
acceptable landing envelope defined as a region in which the relative 
translational and rotational state of the aircraft with respect to the 
pad at touchdown is within some prescribed bounds. Given such a de- 
finition of an acceptable landing, an answer to the second question should 
indicate whether a constant rate-of-descent letdown is adequate. 

Otherwise, a landing guidance system must be developed to assure bringing 
the aircraft to an acceptable touchdown under varying environmental condi- 
tions. 

The purpose of this research effort has been to identify, define and 
analyze factors relevant to the successful landing of VTOL aircraft onboard 
a small aviation ship under adverse environmental conditions. The ultimate 
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objective had been the investigation of alternate landing guidance con- 
cepts for improving touchdown performance over that achievable with a 
constant rate-of-descent letdown. Off-line computer simulation of the 
shipboard landing task was utilized as the principal tool for investigat- 
ing the relative merits of the proposed alternative landing guidance 
concepts . 

This effort represents one aspect of the overall problem of ship- 
board VT0L operations. The ideas and concepts developed in this study 
draw upon the work of previous studies, to some extent. Therefore, key 
sources are referenced throughout this document. 

This report consists of five sections. Section 2 describes the 
various elements that must be integrated into a simulation of the VTOL 
shipboard landing task. Particular emphasis is given to the simulation 
of six degree of freedom (DOF) ship/deck motion using a sum of sinusoids 
representation. The concept of an acceptable landing envelope is defined 
and the touchdown performance of a constant rate-of-descent letdown 
controller is analyzed. Following this, a discussion of alternative 
letdown guidance strategies is presented. Among the options, a go/no go 
letdown decision concept based upon the existence of a ship motion pattern 
of alternating quiescent periods called "lulls" and large motion inter- 
lull intervals called "swells" appears to be the most promising from a 
practical real-time implementation viewpoint. 

Section 3 is devoted to the investigation of lull/swell patterns 
in ship/deck motion. An analytical definition of ship motion lulls is 
presented and utilized in the formulation of a real-time algorithm for 
lull/swell classification. Statistics of lull and swell durations and 
motions are analyzed as a function of the algorithmic parameters. 

Section 4 discusses the improvements in touchdown performance 
obtained with a ship motion pattern based letdown guidance strategy. 


3 



Section 5 summarises the conclusions derived from the study and 
outlines a cockpit simulator test program for verifying the project 
f indings . 



2. SHIPBOARD LANDING GUIDANCE 


In order to develop an off-line computer simulation of the VTOL 
shipboard landing task, it is necessary to have suitable mathematical 
models for the individual elements that characterize the landing 
environment. Therefore, a realistic simulation must include models 
for (a) the ship/deck motion, (b) the airvake turbulence due to the 
ship’s superstructure, (c) the aircraft dynamics including the character- 
istics of the associated flight control system, (d) the shipboard and 
airborne navigation system, and finally, (e) the landing guidance 
system. However, the principal objective of this study has been the 
analysis of VTOI, landing guidance system requirements to ensure 
acceptable impact conditions at touchdown onboard a moving ship under 
rough sea conditions. The primary factor which makes shipboard lancing 
sc difficult is the existence of large and apparently unpredictable 
excursions in landing pad position (primarily heave) and orientation 
(mostly rcll and pitch deviations) caused by the wave disturbances in 
a high sea state {>_ 5> environment. Therefore, greatest emphasis was 
placed during this study on realistic ship motion simulation. Ship 
airwake turbulence and deck proximity effects are not considered. 
Furthermore, it is assumed that the airborne equipment has access to 
exact: (a) ship attitude, (b) aircraft attitude, and (c) aircraft and 
landing pad position coordinates in a local level reference frame 
centered at the average position of the ship's center of gravity or 
bullseye. The aircraft control system is assumed to provide first or 
second order step response characteristics in the angular and transla- 
tional degrees of freedom. Details of the system bandwidth and time 
constant parameters are discussed in section 2. A where alternate landing 
guidance and controller concepts are considered. 


2.1. Ship Motion Simulation 

The NAVAIRENGCEN Ship Motion Computer Program [1] was used in this 
study to generate ship motion time history data. The program is based 
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on a technique for modeling ship motions that was introduced by St. Denis 
and Pierson in their classic paper "on the Motions of Ships in Confused 
Seas" [2]. The approach consists of representing the ship motion in each 
of the six degrees of freedom (roll, pitch, yaw, surge, sway, heave) at 
the ship's renter of gravity ( CG) by a sum of sinusoids. Selection of 
the number of sinusoids and their amplitudes, frequencies and relative 
phase angles is based upon having a dynamical ship model, as shown in 
Figure 2, for computing the ship motion response to wave disturbance 
inputs. Ship motion response h(t) is related to the wave input C(t) throucl 
-he ship transfer function (also referred to as the regular wave transfer 

function in Ref. 1) H (jw) . 

n 


Seaway 

Forcing 

Function 


Ship Transfer 
Function 


n(t) Ship 

— — ► Motion 

S (co) Response 



H r (jw) 

III 

5(t) 

S ; (u) 

“J $_,(“>) 
|H (jco)|e 1 


Figure 2. Ship Motion Response Computation 


The seaway or wave disturbance ;(t,> is defined to be a stationary 
ergodic random process with power spectrum S^(w) , where w is the frequency 
(rad/s) of the wave inputs. Note that the definition of spectral density 
in wave theory is slightly different than the standard form used in 
linear systems theory. Thus, spectral density in wave theory is defined as: 


s x (.) 


A 

A 

m 



y T > e 


-jU)T 


dT 


0 


cc _> 0 

OJ < 0 


(1) 


Similarly, autocorrelation is defined as 


R x (t) - J Re[S x (aOe jWT ]dw 

0 


(2) 
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Thus, mean square value of x(t) is 




S O) dw . 
x 
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( 3 ) 


For the purpose of this study, 
used to represent the wave spectra. 


the Bretschneider spectrum [1] is 
Thus , 


_ , . A A 
s c (^) » ~r exp 


('i) 


(4) 


where 


483.5 

T 4 
0 

1944. 5 


’1/3 


'1/3 


significant wave height (average double amplitude of 
the one third highest waves) 


and 


A 0 


modal period. 


Figure 3 shows the normalized Bretschneider spectrum for . rious 

values of Tq. Note that the modal period car. be used to adjust the 

peak magnitude frequency oj according to the relation 

max 


w 

max 


6.28 


rad/s 


( 3 ) 


This relationship may be used to tune the ship motion spectrum in any 
degree of freedom. 


The ship motion response spectrum in any degree of freedom is 
computed as 


S (w) - lH (ju) 2 S r (u) 

m n 


( 6 ) 
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Figure 3. Normalized Bretschneider Spectrum for Varying 
Modal Period Tq. [From Ref. 1] . 


In order to compute the response spectrum in the encounter frequency 
(i~ e ) domain, it is necessary to use the mapping 


OJ 


a 


f (60) = uJ - — V COS u 

g s s 


(7) 


whe re 


V = ship speed 


p * ship heading with respect to the predominant 
wave direction 

2 

g = acceleration due to gravity (32.2 ft/s ). 


Thus , 


W 


V") 


oOu 

e 

oCG 


■d U) = f ( GJ f 


V“> 


1 - ~ V cos p 
L g s H s 


-* CJ * 


f" 1 (W o ) 

e 


( 8 ) 
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where 
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w « f~% e ) 


V 


-1 +^ 1 + 


/ 4^ V cos u 

e s s 


2 V cos U 
s s 


Figure 4 illustrates the steps involved in computing S (u> ). 
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2 
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Figure 4. Ship Response Spectrum Computation 


(9) 


In wave theory, the square of the ship transfer function magnitude is 
called the Response Amplitude Operator (RAO). Thus, 

RAO^co) = | H^C jo>) j 2 (10) 

RAO data for each of the six degrees of freedom are usually generated 
by measuring ship response in towing tank experiments or derived 
computationally by exercising the NSRDC Ship Motion and Sea Load Computer 
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.• - ■ ■■ - ]•* :/ 

r r v.‘ , . * ^ - -is i 

Program [3]. It should be noted that RAO data depend upon the type of 

ship, the ship’s speed (V ) and heading (y ), and the degree of freedom. 

s s 

The sum of sinusoids model for describing ship motion time history 
for any degree of freedom is expressed as: 


n(t) = n t cos (w ei t - 

i=l 


( 11 ) 


where N = number of sinusoids 

ft) ^ = location of the i th frequency component 


$. ■ - H (jcu. ) 

i q J i 


co. = f _1 (co .) 
l ei 


<j> (to . ) 
n i 


= phase lag of the ship transfer function at 
the i th frequency 

H i = amplitude of i _th frequency component 

£. * random phase of i _th frequency component 
1 ~ [0 - 2tt] 


• | H (ju) 

i 1 n_ 


Note that: jco) = | ( j co> j e 


= VRAD (w) e 


- j $ (co) 

n 


The N frequencies ft) are equispaced 5co e apart over the encounter 

frequency range [co . , co ] where the ship motion spectrum S (co ) 

J ° e min e max v r H e 

exceeds 4% of its maximum value. The corresponding N amplitudes ri^ are 
chosen such that the power of the sinusoid at ft) approximates the area 
under the power spectrum S^co^) over the frequency bin 


6<o 


co 


ei 


6 co 1 

CO + 
ei 2 


Thus , 


’•i A 

mr * S (ft) J 6ft) 

2 n ei e 


( 12 ) 
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or 


v/TyS-T^ET 


(13) 


Note, however, that the NAEC Ship Motion Computer Program [1] used 
in this study, performs spectral deconposition into a sum of sinusoids 
in the wave frequency (w) domain. In other words, K equispaced frequencies 
u). are selected in the m domain with computed as 

n ± = / 2 • S n (0O • 6(0 (14) 

where 6u> is the constant frequency spacing. Note that, calculated using 

Eq. (14) is identical to that obtained from Eq. (13) with 6w ei substituted 

for ooj e . These amplitudes{n^T- are assigned to frequencies in the 

encounter frequency domain, calculated using the frequency mapping 

defined in Eq. (7) between u and . The incremental frequency bins 

6u . are related to bv the nonlinear Jacobian transformation, 
ei 




6<j 


(15) 


Consequently, the frequencies in the domain are not uniformly spaced. 

In fact for the case (V = 25 kt , p = 120°) considered in this study, 

S S 

£ 

the Jacobian - — over the frequency range of interest is slightly 
greater than one. The result is a slow rise in bin size Co.' e ^ as i 
increases from 1 to N. However, for the purpose of this study, this 
effect was not considered to be serious enough to warrant a modification 
of the NAEC program. 


2.1.1 Generation of Simulated Ship Motion Data 

The NAEC Ship Motion Computer Program [1] was used to generate 
simulated ship motion time histories. However, in order to use the 
program, it is necessary to specify the following three sets of 
parameters : 

a) ship type, speed V g and heading y ; 
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b) sea spectrum parameters - namely, significant wave height 

^1/3 modal period T Q ; and 

c) ship dynamics parameters - namely the RAO^w) and phase 
angles (^(w) over the frequency range of interest, for 
each degree of freedom rj. 

The DD963 Spruance class destroyer data base provided by the NAEC program 
was used to generate the motion time histories. The program data base 
covers an operating range of ship speeds V g from 5 to 25 kt in 5 kt 
increments, and ship headings from 0 (following seas) to 180 (head seas) 
in 15 degree increments. Figure 5 shows the parameters characterizing 
the ship motion relative to the incident wave direction. For the 
purpose of this study an operating condition corresponding to V g = 25 kt 
and ti g = 120° was selected because it results in large deck motions in 
all degrees of freedom (in particular, heave and pitch). 

The significant wave height of 32 ft (9.75 m) was chosen to 

correspond to a high sea state condition. The modal period Tq was 
selected by a tuning process as illustrated in Figure 6. The objective 
of the tuning process was to select a value for Tq that results in a heave 
deck motion spectrum having narrow band characteristics that are typical 
of actual ship motion. Figures 6a and 6b show the peak magnitude 
normalized sea spectrum, deck heave RAO, and deck heave spectrum for 
Tq = 13.5s and 8.48s, respectively. The value of Tq = 8.48s tunes the 
peak of the sea spectrum with the peak of the heave RAO, giving rise to 
a narrow band heave motion spectrum. 



Figure 5. Parameters Characterizing Ship Motion Relative to 

Incident Wave Direction. 
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O Sea Spectrum S (m ) 
£ e 

□ Heave RAO (jj ) 
e 

A Deck Heave Spectrum 
^Heave^ a> e'* 



Figure 7 shows the simulated DD963 ship motion time histories using 
a sum of seventy sinusoids over a 1200s interval for the six deck motion 
degrees of freedom. The motion is described in a local level reference 
frame centered at the average position of the bullseye, with x axis 
astern, y axis starboard and z axis up. The deck motion quantities shown 
are : 


x - longitudinal 
y - lateral 
z - heave or vertical 
<t> - roll 

8 - pitch 

\p - yaw 


( 16 ) 


The ordinate scales for dimensionally equivalent translational and 
rotational quantities are identical to facilitate visual comparison of 
their relative magnitudes and temporal characteristics. Examination of 
these time histories reveals that the motions display temporal character- 
istics that are typical of narrow-band random processes. Specifically, 
the time histories show an alternating sequence of quiescent and large 
ship motions, referred to in the literature as lulls and inter- lulls 
(or swells) [ 4 J - Furthermore, it appears that the lulls and inter-lull 
intervals are synchronized across the six degrees of freedom. 


However, previous work by Fortenbaugh [5] concluded that a sum of 
six sinusoids (N = 6) should provide an excellent approximation for 
representing actual ship motions. This conclusion was based on showing 
that root-mean-square (rms) position, rate, and acceleration for each 
of the six degrees of freedom differ by less than 5% for the six or 
thirty- two sinusoids approximations. This agreement of rms values 
does not guarantee equivalence of the underlying distributions, 
correlation functions, and simulated time histories of the motion in each 

degree of freedom. This is apparent on a cursory examination of the time 

• • • 

histories for [<{>, <£, 0, 0, z, z] for the six and seventy sum of sinusoids 
models, shown in Figures 8 and 9, respectively. Most importantly, the 
time histories (Fig. 8) for the sum of six sinusoids do not display the 
alternating lull and swell periods so obvious in the seventy sinusoids 
model (Fig. 9 and Fig. 7). Furthermore, the six sinusoids data has a 
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Figure 7. Simulated Ship Motion Time Histories (70 Sinusoids): 
V g - 25 kt, p g - 120°, c 1/3 - 32 ft, T Q - 8.48 s. 
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Figure 7. (Cont.) 
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Figure 8. Simulated Ship Motion Time Histories (6 Sinusoids) 
V g - 25 kt, u 9 - 120°, ; 1/3 « 32 ft, T Q = 8.48 s. 
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Figure 9. Simulated Ship Motion Time Histories (70 Sinusoids) 
V - 25 kt, y - 120°, C,., - 32 ft, T rt - 8.A8 s. 




definite repetitive pattern with a fundamental period of approximately 
340 seconds. Differences in the underlying distributions for the two 
cases are also noticeable in the histograms for deck heave (z) and 
vertical velocity ( 2 ) shown in Figure 10. The sum of seventy sinusoids 
distributions are markedly more Gaussian than the ones for six sinusoids, 
the most noticeable discrepancy being in the six sinusoids histogram 
for z. The whole issue of developing minimal complexity models (including, 
but not limited to sum of sunusoids or state-space models) for real time 
ship motion simulation is an interesting area for more systematic research 
but. was considered to be outside the scope of this study effort. 

Henceforth, ship motion data generated by the NAEC Ship Motion 
Computer Program using a sum of seventy sinusoids representation (see 
sample time histories in Figure 7) will be used as a baseline for 
investigating shipboard landing guidance system c"' r >^' a Dts. 


2.2 Landing Envelope Definition 


Landing envelope is a term used to define the acceptable set of 

conditions at touchdown. Typically, these requirements are stipulated 

in terms of the aircraft state x , the landing pad state x_ , and the 

A V A 

relative state of the aircraft with respect to the pad x^j p(* x A *x p ) 
at touchdown. Note that the state vector is assumed to be composed of 
the three position (x,y,z) and three attitude (ty, 6,<$>) components, and 
their first and second derivatives . Thus, given the constraints on 


a. > a 

— -p 


anu x 


-A/ T 


as 


and 


—A 




£ S, 


e S, 


-A/P G S A/P 


( 17 ) 

(18) 
(19) 


where "e " means "belongs to" or "is a member of" and "S" denotes a "set" 
in state-space coordinates. 


20 



rRCOUL'NCY PRCOliCNCT 


ORIGINAL PAGE 19 
OF POOR QUALITY 


70 SINE WRVES 70 SINE WRVES 



2 w 2 r*/% 


6 SINE NAVES 6 SINE WAVES 



2 m 2 / £ 


Figure 10. Probability Density Functions for z and z for Sum of 

70 and 6 Sinusoids. 


21 



the landing envelope E is defined as the set of aircraft and landing pad 
states satisfying Eq. ( 17) — ( 19) . Thus, 


E 


£ 

s 



x . e S . and e S-, 
—A A — P P 


a„d x A/p 



( 20 ) 


A review of previous work on the shipboard landing problem indicates 
that there is no consistent definition of a landing envelope in the 
sense defined by Eq. (20). However, it is generally agreed that, as a 
minimum, the relative attitude and translational velocity of the air- 
craft with respect to the landing pad at touchdown should be within 
reasonable bounds. For example, at the time of touchdown, conditions 
such as the following may be stipulated: 


A/P 

< 

$ 

max 

A/P 

< 

e 

max 

A/P 

< 

max 

'A/P 

< 

• 

X 

max 

A/P 

< 

• 

y 

max 

A/P 

< 

• 

z 

max 


( 21 ) 


Unanimity on selection of the numerical values for these limits is 
not possible because of their dependence on ship type, aircraft landing 
gear characteristics and other operational factors. Additional limits 
on landing pad state at touchdown are usually imposed because of 
limitations in the aircraft control power (thrust /weight ratio) , and the 
speed cf response (i.e. bandwidth) of the attitude and vertical transla- 
tional flight control systems. These aircraft bandwidth limitations 
preclude the design of a ship motion following landing control system 
during constant speed letdown. As a result, previous studies [6-8] nave 
proposed the following requirements bn the landing pad motion at touchdown 
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a) 


b) 


c) 


The pad vertical position with respect to its average value 
shall be greater than zero. 


i .e. : 


z > 0 


( 22 ) 


The vertical velocity of the pad shall be positive (moving 
upward) but shall be less than some prescribed value (e.g. 
.6 m/s or 2 ft/s) 


i .e . : 


0 < z < z 


max 


(23) 


The pitch and roll attitudes and rates shall be less than 
some prescribed values 


i.e. : 


'e < e 

1 — max 


! e I < el 


max 


(24) 


1<J>| 1 4>j 
ill < +1 


max 


max 


Equations (22)-(24) represent one plausible set Sp (see Eq. 18) of acceptable 
landing pad states at touchdown, and provide an example for elucidating the 
underlying concept. These three constraints reflect a conservative landing 
guidance philosophy where the landing window is limited to periods during 
which the ship’s landing pad motion is small and at its peak or crest values. 
There restrictions are imposed to minimize the probability of large heave 
velocity at touchdown thereby reducing the chances for a hard landing. 


The basic objective of landing guidance is to bring the aircraft to 
a safe touchdown onboard the deck while staying within its control system 
limitations. A careful review of previous work on the subject indicates 
little agreement on the parameters defining an acceptable landing envelope. 
As a result, a normative approach is used in this study towards evaluating 
alternative landing guidance systems concepts, whereby improvements in 
touchdown performance obtained with alternate designs are normalized with 
respect to 1 ae performance of a simple constant-rate-of-descent letdown 
controller. 
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2.3 Constant Rate-of-Descent (CROD) Letdovm 


ORIGINAL PAGE IS 
OF POOR QUALITY 


The simplest form of letdown procedure is to use no active guidance 
by allowing constant rate-of-descent to a touchdown at any time after the 
aircraft has achieved a stable hover position over the ships landing pad. 
The touchdown performance achieved with the CROD letdown profile may, 
therefore, be examined to establish a baseline or standard of performance 
against which alternative letdown guidance schemes may be measured. 


Figure 11 illustrates the CROD letdown scenario. Throughout the 
remainder of this report, the stationkeeping or hover point is assumed 
to be at a constant altitude cf 20 ft over the average position of the 
bullseye. This height was chosen to allow a safe clearance above the 
maximum deck height (z 15 ft) obtained under the simulated landing con- 


ditions (i.e., V * 25 kt, y 
s s 


120°, 5 1/3 - 32 ft, T Q = 8.48 s) 


CONSTANT RATE Of DESCENT I2ft/s) TOUCHDOWN 

"•°l T r - Hover Altitude T 







m/s 



Note that for a 2 ft/s letdown rate, touchdowns occur only when pad 
vertical velocity is greater than -2 ft/s (i.e. above the cross-hatched 
region in z) . Since the principal factor limiting shipboard landing 
is the relative vertical impact velocity between the aircraft landing 
gear and the ship's deck at touchdown, only the landing pad vertical 
motion (z) and vertical velocity (z) are shown over a 100 second 
segment. Two CROD (2 ft/s) letdown profiles are illustrated. The 
figure shows that profiles 1 and 2 land during different periods of 
ship motion. Letdown profile 1 lands (T^ * 578s) when ship motion is 
in a relatively quiescent period while profile 2 touches down (T^ “ 608 s) 
during 'arge amplitude period. As a result, the touchdown vertical 
impact velocity is 1 ft/s for profile 1 and 8 ft/s for profile 2. The 
CROD letdown is performed with the aircraft attitude at trim value 
(assumed to be zero attitude throughout this report i .e. ip^ = 9^ = = 0°). 

Therefore, touchdown performance in the remainder of the translational 
(i.e. x, x, x, y, y, y) and rotational ($,'<£, <P, 9, 9, 9, ip, $) 
variables may be expressed directly in terms of the landing pad conditions 
at touchdown. 

Touchdown performance can be evaluated by examining the histograms 
(or probability density functions - PDF's), cumulative distribution 
functions (i.e. F(x) = Pr [X <_ x ]) or survival functions (i.e.t^x) = Pr 
[X > x] = 1 - Pr [X <_ x] = 1 - F(x)) for the individual variables. A 
Monte Carlo simulation approach was used in this study to compute the 
touchdown statistics. The simulation allowed for a letdown profile to 
be initiated every 1 second. Touchdown data over a one hour period 
(approximately 3600 touchdowns) we re used to compute the histograms and 
distributions of the individual state variables. An examination of the 

ship motion time histories in Figure 7 indicates that out of the eighteen 

• • ta • 

variables shown, only four, namely z, y, y and <p , may be limiting factors 
in achieving an acceptable touchdown performance. Therefore, discussion 
of the results will be limited to these four quantities. Touchdown 
performance will be characterized in terms of survival functions for the 

• • ft • 

vertical impact velocity z R (H z - z^) , lateral impact velocity y, lateral 
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impact acceleration y (proportional to side force on landing gear) , and 

• • 
roll angular velocity 4>. Plots for survival function are given for z , 

alone. However, survxval function results in a tabular format are 

provided for each of the four variables. Note that the survival 

function gives the probability of a given variable exceeding a prescribed 

threshold. This function, therefore, provides a direct measure of the 

touchdown performance in concrete and tangible terms. 

The survival function, t^(z ) , for the 2 ft/s CROD letdown guidance 

K a # 

scheme is shown in Figure 12* Corresponding results for z , ]yj, }y| and 
• ^ 

|$| are given in a tabular form in Table 1. Note that an absolute value 

on z c is not necessary since z_ >0, by definition. Note that 
K K 

^(| y|) ■ Probability that absolute lateral impact 

velocity at touchdown is greater than |y| ' ' 

The same definition applies to ^(|y|) and &(\<t>\) . Touchdown performance 
results for the CROD letdown profile summarized in Table 1 shall be used 
throughout the remainder of this report for evaluating, on a normative 
basis, the improvements over CROD guidance provided by alternative 
letdown schemes. 


2.4. Alternative Letdown Guidance Concepts 

The performance of a CROD letdown scheme is summarized by the results 
in Table 1. The purpose of this section is to discuss alternative letdown 
guidance concepts which give improved touchdown performance by decreasing 
the probability of touchdowns with large translational and rotational 
velocities and accelerations at impact. 

Two categories of letdown guidance laws may be formulated: 

1) Designs based upon the assumption that ship motion response 

time histories display a pattern of alternating low amplitude 
quiescent periods called "lulls" followed by large motion 
inter- lull segments called "swells", and 
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TABLE 1. 


ORIGINAL PAGE !9 
OF POOR QUALITY 


CROD Letdown Touchdown Performance 


N • 
TO 

i 

• 

y 

y 

\ ^R> 

(ft/s) x 100% 

1*1 t^(|$|) 

(deg/s) x 100% 

i y ! «^(!yj) 

(ft/s) X 100% 

I y | [ y ! ) 

ft/s/s x 100% 

6 44 

8 22 

10 9 

12 3 

2 44 

3 23 

4 10 

5 3 

2 54 

4 19 

5 9 

7 2 

2 43 

4 12 

5 5 

7 1 


2) Designs that are independent of the presence or absence 

of a ship motion pattern structure (i.e. lulls and swells). 

Previous studies [6-8] of the shipboard landing problem have ignored 
the existence of ship motion lulls (i.e. design option #1) and have instead 
proposed two approaches to letdown guidance: A ship motion following 

controller and a ship motion forecasting controller. The CROD letdown does 
not allow the aircraft to respond to instantaneous ship motion, and in 
particular, to its relative attitude and vertical velocity with respect to 
the moving landing pad. One approach to reducing the relative motion at 
touchdown is to design a flight control system or landing controller that 
forces the aircraft to track or follow the instantaneous landing pad motion 
while maintaining a CROD letdown until touchdown. The basic idea is 
illustrated in Figure 13 which shows the application of ship motion following 
letdown guidance to the vertical motion degree of freedom. 
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Figure 13. Ship Motion Following Letdown Guidance for Vertical Motion. 


Typical response of this design is illustrated in Figure 14. The total 
vertical velocity command z is the sum of the nominal CROD command 



Figure 14. Response of Vertical Ship Motion Following Controller. 
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• _ • 

^ROD^~ ” 2 ft/s) and a feedback control z . Thus, 

Z T " Z CR0D + z c (26) 

The command z is provided by a landing controller whose goal is to 
minimize the instantaneous position and velocity errors in following 
vertical ship motion z. The ship motion following landing controller may 

be designed using classical and/or modern control theory methods; the 
choice depending upon the nature and quality of the apriori information 
(i.e. mathematical models for the aircraft, landing pad motion and exter- 
nal gust disturbances) . Such a letdown controller can be used only if 
the aircraft can follow the landing pad motion with reasonable fidelity. 
This would require a flight control system with a sufficiently high 
bandwidth (at least three times the dominant frequency in the given ship 
degree of freedom ) , servo actuator velocity and acceleration authority 
limits for successful implementation. However, existing aircraft do not 
have such characteristics and hence are not capable of following the 
relatively high frequency (0.3- 2.0 rad/s) deck motions. As a result, 
ship motion following letdown guidance is not considered to be a viable 
technique for shipboard landing [8]. 


An alternative approach is to use some kind of ship motion fore- 
casting algorithm to determine the predicted time periods or landing 
windows during which aircraft landings would be acceptable. A suitable 
letdown guidance scheme (e.g. CROD or LQG control with terminal constraints 
[7]) may be used to bring the aircraft to an acceptable touchdown during 
the designated landing window. The specific criteria for designating a 
future time period as a landing window would depend upon the touchdown 
requirements stipulated by the acceptable landing envelope. However, 
developing a reliable ship motion forecasting algorithm is a non trivial 
problem because the efficacy of forecasting methods is predicated upon the 
veracity of analytical models (e.g. finite order state space or auto- 
regressive-moving average (ARHA) formulations) used to represent the ship 
motion spectra (auto and cross spectra). Although the subject of time 


3C 



series forecasting has received considerable attention over the years [9], 
and continues to be an area of ongoing research, the usefulness of the 
available model-based techniques for extrapolation or forecasting of real 
(as opposed to computer generated) time series is, at best, questionable. 
Models for ship motion response need to be developed and validated using 
systematic structure determination and parameter identification techniques 
with real data obtained from sea trials. However, it must be emphasized 
that even with exact models for given time series data, the time window 
over which forecasts can be made with a high level of confidence is 
finite and limited due to the stochastic nature of the underlying processes. 

A ship motion pattern directed letdown guidance scheme was developed 
in this study. The approach exploits the fundamental temporal character- 
istic of alternating lull/swell patterns so evident in actual ship motion 
data and the sum of seventy sinusoids representation of Figure 7. A real 
time algorithm for lull/swell classification based upon ship motion 
pattern features is described in the next section. The output of the 
classification algorithm is used to command a go/no go signal to indicate 
the initiation and termination of a landing window. Performance improve- 
ments obtained with such a go/no go pattern based letdown guidance scheme 
are discussed in section 4. 
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3. SHIP/LANDING PAD MOTION PATTERN ANALYSIS 


Ship motion response time histories as shown in Figure 7» can be 
classified into one of two categories or modes: (1) a "lull" mode 
representing time intervals where overall ship motion may be described 
as relatively low, and (2) a "swell" mode where ship motion is relatively 
high. Since the lulls and swells for all the six degrees of freedom 
appear to be synchronized, any one of the eighteen variables in Figure 7 
may be used for lull/swell categorization. 1° this effort, the landing pad 
vertical motion z(t) was used to develop a real time lull/swell classifica- 
tion algorithm. 


The concept of a lull is best described by reference to Figure 15 
where the top two plots (a) and (b) show sample landing pad vertical 
position z(t), and velocity z(t), over a one hundred second interval 
(identical to z and z in Figure 11). The lull/swell classification 
algorithm is characterized by an indicator function I(z) which maps the 

• it 

ship motion z(t) (including, if necessary, higher derivatives z and z) 
into a binary output sequence of zeros and ones, where I(z) » 0 corres- 
ponds to a lull period and I(z) * 1 describes an inter-lull interval or 

• „ 

swell. An examination of the vertical motion time histories (i.e., z, z, z) 
in Figure 15 (and Figure 7 for a 1200 second period) shows that the most 
obvious way to dichotomize the motion into lulls and swells is by 
defining an indicator function in terms of a threshold amplitude z^, such 


that 


( 0 (Lull) : z(t) <_ z T 
( 1 (Swell) : z ( t ) > z T 


(27) 


However, the lull/swell switch plots 1^ that are obtained using this 

definition display a high switching frequency (0 to 1 and vice versa) 
making such a criterion unacceptable from a practical viewpoint. The 
reason for this behavior is the existence of a relatively large number 
of single peaks in z(t) exceeding (below) the threshold z^. that occur during 
otherwise long periods, where z(t) <_ Zj ( z(t) > z^ ). The net result is 
a large number of lulls and swells of extremely small duration (<_ 3 s) 
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chat are imbedded in relatively longer duration swell end lull intervals. 

A more reasonable approach is to define the initiation and 

termination of a lull in terms of parameters that reflect the motion 

response characteristics over a finite past or memory interval. Thus, 

initiation of a lull may be defined by a parameter "NSTART" as follows. 

A lull is said to begin (i.e.: I switches from 1 r.o 0) if NSTART 

z 

consecutive positive peaks of z(t) fall below a prescribed threshold 
value Zj,. Similarly, termination of a lull (or onset of a swell) may 
be defined by a parameter "NSTOP" as follows. A lull is said f.o have 
terminated (i.e. 1^ switches from 0 to 1) when NSTOP consecutive 
positive peaks of z(t) exceed a prescribed threshold value z^. Figure 
15(c) shows a plot of the indicator function I z for the lull/swell 
periods for nominal parameter values of NSTART ■ 2, NSTOP ■ 2 and z j ■ 

5.04 ft. The threshold value of Zj was chosen as the mean value of the 
positive peak cmplitude envelope of z(t). Over the 100 s segment of 
data shown in Figure 15, the interval (602.2 - 644. 4)s is identified as 
a swell period. In other words, the. ongoing lull from t * 555 s is 
terminated at t • 602.2 s whan two consecutive positive peaks of z(t) 

(i.e. NSTOP » 2) have exceeded the threshold of Zj ■ 5.04 ft (see Figure 
15(a,c)). Similarly, a lull is initiated at t * 644.4 s when two con- 
secutive positive peaks of z(t) (i.e. NSTART * 2) have stayed under the 
threshold of z^ ■ 5.04 ft (see Figure 15(a,c)). 

Note that NSTOP « 2 results in true swells (i.e. swells which 
usually last for >3*4 cyclec or approximately >_ 20 seconds) being 
detected one cycle after their actual onset. Selecting NSTOP » 1 would 
apparently solve this problem by detecting true swells one cycle earlier 
than with NSTOP ■ 2. However, this would create too many additional swells 
of 6 - 10 seconds duration due to the occurrence of isolated single peaks 
of z(t) exceeding the prescribed threshold Zj. The termination of lulls 
can be made one cycle earlier, when the first peak of z(t) exceeds z^., 
if at that time one can confidently forecast that the next successive 
positive peak shall also exceed the prescribed threshold. This is illus- 
trated in Figure 15(a,c), where the one cycle ahead forecasting based lull 
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termination occurs 5.4 seconds earlier it t ■ 596.8 8. Unfortunately, a 
reliable on-line one cycle ahead forecasting method is needed for practical 
implementation of the approach. 


An alternative method for early detection of the onset of a swell 

e 

based upon the z versus z phase plane pattern characteristics at lull/ 
swell transitions was developed during this study. Phase plane plots 
of s(t) versus z(t) for several lull/swell transition segments (as 
defined by NSTOF * 2 and z^, * 5.04 ft) were analyzed to determine if 
a pattern or a sequence of events in the z and z time histories could 
be identified and .sed as a precursor for lull termination prior to the 
crossing of the threshold z^. by the second positive peak of z(t). 

Based on thif> analysis, a phase plane based criterion for lull termina- 
tion that performs nearly as well as the one peak ahead forecasting 
based method was developed. According to this critei*un, a lull is 

terminated if successive peak magnitudes of vertical position z(t) and 
• • 

vertical velocity z(t) (or z followed by z) exceed their prescribed 

thresholds. Thus, I *, the indicator function for this phase plane 

z ,z 

approach switches from zero to one according to Eq. (28). 


. £ 


Z.z 


l 2 Peak (t) 


> z T and 


"Peak 


( t ) { > z„ 


or 


2 Peak (t) 


> Zj and 


"Peak 


(t )i > z. 


(28) 


« ^ ^ • 
wher^s 2 Poa’, (t > (z Peak (t )} is the first peak in z(t) (z(t)) following 

ih? event |z pe _^(t) [ > z T (Up eak (t) ! > z p ) . Figure 15(d) shows a plot 

of the resulting indicator function I •. Notice that the end of the 

z.z 

lull is detected at t * 596.4 s. This compares favorably to the per- 
formance of the one peak ahead forecast based threshold test shown in 
Figure 15(c). This is further apparent on comparison of the indicator 

functions I and I , • over a period of 1200 seconds as shown in 
z z z 

Figure 16. 


The effects of varying NS TART (1,2 or 3) on I • and hence 
luil/svell classification are shown in Figure 15 (e,d,f), respectively. 
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Figure 16. Comparison of I and I , • over a 1200 Second Period 





As expected, increasing NSTART values have an impact on lull/swell 
duration statistics. Figure 17 shows the survivor functions ^(TL) and 
^(TI) for the lull duration and inter-lull (swell) interval, respec- 
tively, for NSTART =1,2 and 3. These plots show that the principal 
effect of increasing NSTART is on the swell interval statistics, and 
most significantly on the likelihood of swell periods of two minutes 
and longer. A go/no go guidance command based upon a lull/swell classi- 
fication algorithm with such parameters (NSTART = 3) could result in a 
no-go command lasting for two minutes or longer, which would be un- 
acceptable because of constraints imposed by pilot workload, fuel 
reserves, and other operational considerations. Therefore, either 
NSTART = 1 or 2 may be used to define the initiation of a lull. A value 
of NSTART = 2 was selected as nominal during this study. 

• • • 

Figure 18 shows the histograms of z, y, y and ? during lull and swell 
periods corresponding to the I 2 ^ phase plane indicator function shown in 

Figure 16(d). These histograms show that the standard deviations 

• • • 

of the variables z, y, y and 4> have significantly lower values during lull 
periods in comparison to their values during swell intervals. Therefore, 
improvements in the touchdown performance should be possible if shipboard 
landings can be restricted to occur during lull periods alone. 

The lull and swell duration survival functions^(TL) ande^TI) for the 
three options of lull termination shown in Figure I5(c,d), namely a) z thres- 
hold, b) z threshold with one peak ahead forecasting, and c) z,z phase plane, 
are plotted in Figure 19. The plots show that the options (b) and (c) above 
are nearly the same. However, they give slight lower (higher) values of the 
lull (swell) duration survival functions than option (a). 

In summary, a generic structure for an on-line/swell classification 
algorithm has been developed. The selection of the threshold values z^, and 
z^, and parameter NSTART define the lull/swell classification rule and the 
resulting temporal and statistical characteristics of the lull and swell 
patterns. 
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Figure 17 , Survivor Functions ^(TL ,) and t^(TI) for Lull Duration 

TL and Swell Interval TI. 
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Figure 18. Lull/Swell Amplitude Histograms for z and z with 

I z • and z T = 5.04 ft, z^ ■ 7.5 ft/s, NSTART = 2. 
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Figure 19. Lull and Swell Duration Histograms for Three 

Lull Termination Options. 
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The next section is devoted to an evaluation of alternative pattern 
based letdown guidance schemes. Improvements in touchdown performance 
obtained using a lull directed adaptive letdown guidance strategy are 
presented . 



4. SHIP MOTION PATTERN DIRECTED LETDOWN GUIDANCE 


The purpose of this section is to investigate the potential improve- 
ment in touchdown performance that can be achieved with a lull/swell 
pattern directed letdown guidance strategy. Three different lull/swell 
classification schemes described in the previous section are considered. 
They are: 

Scheme (a) lull/swell classification based upon the use of a heave 

amplitude threshold (z^,) driven indicator function 

I with z« 5.04 ft, NSTART = 2 and NSTOP = 2 
z i 

(see Fig. 16, dark switch plot in l ^) ; 

Scheme (b) lull/swell classification based upon the use of indicator 
function I as in (a) above but with the additional capa- 
bility of one peak ahead forecasting or prediction (see 
Fig. 16, dotted switch plot in I ); and 

Scheme (c) lull/swell classification based upon the use of a deck 
vertical motion phase plane threshold (z^,, 7. m ) driven 

_ x 

indicator function - I , • with ?_ = 5.04 ft, z_ = 7.5 ft/s, 

z z 1 x 

and NSTART - 2. 

A CROD (-2 ft/s) letdown profile may be initiated at any time during 
a lull period. Three options are available at the termination of a lull 
with respect to the future trajectory of ongoing letdown profiles, as 
follows : 

Option (1): Continue ongoing CROD letdown profile until termination at 

touchdown . 

Option (2): Arrest the aircraft rate-of-descent and hold the altitude 

of the aircraft at lull termination until touchdown. 

Option (3): Abort the aircraft rate-of-descent, initiate maximum rate- 

of-climb, and return to the prescribed hover altitude. 
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Eight alternative letdown guidance strategies defined in Table 2 were 
simulated and analyzed. Strategy #0 represents the nominal CROD letdown 
guidance profile discussed earlier in Section 2.3. An off line computer 
simulation of the landing task was conducted to generate the touchdown 
performance statistics. Letdown profiles from the hover point until 
touchdown were initiated every one (1) second when permitted to do so by 
the particular guidance strategy. Touchdown performance is presented in 
terms of the survival functions for z , j | and |y| . The results are pre- 
sented in two formats - plots and tables. Figure 20 shows the survival 

function ) for letdown guidance strategies 0 through 3, respectively. 

K 

Similarly, Figure 21 shows <^(z ) for letdown guidance strategies ?/0 , #4, 

#5, #6, and #7. Survival functions for z^, |y| and j$| are tabulated in 
Tables 3-5, respectively. 

Figure 20 shows that there is no significant difference in the survival 
functions <^\z R ) for the four strategies 0 thru 3. In other words, no improve- 
ment is obtained by using a letdown guidance strategy wherein the decision 
to initiate a letdown ("go" decision) is restricted to lull periods alone. 

This is because the problem lies at the termination of a lull with the 
letdown profiles already initiated in the past ten (10) seconds and which 
have not yet landed. The result is that these profiles terminate in touch- 
downs during swell ("no go" decision) intervals, when deck vertical velocity 
z can be high. Thus, for guidance strategies 1 through 3, although the actual 
number of touchdowns with large relative impact velocities (e.g., z > 12 ft/s) 
is lower, the percentage of these touchdowns with respect to the total number 
of touchdowns is not significantly different from that for the nominal CROD 
letdown strategy #0. 

Letdown profile options (2) and (3) provide a way for handling these 
ongoing landings at the instant of lull to swell transition. In option (2), 
the number of touchdowns during a swell interval, as well as the relative 
vertical impact velocities at touchdown, is reduced by switching the let- 
down rate-of-descent from -2 ft/s to zero at lull termination. Thereafter, 
the aircraft is maintained in an altitude hold mode until touchdown during 
a swell or until the initiation of the next lull period ("go" decision) when 
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TABLE 2 

Alternative Letdown Guidance Strategies 


Letdown 
Guidance 
Strategy if 

Lull/Swell 

Classification 

i 

Scheme 

Letdown Profile 
Option at Lull 
Termination 

0 

None 

N/A 

CROD Letdown 

1 

(a) 

(1) 

2 

(b) 

(1) 

3 

(c) 

(1) 

4 

(a) 



(2) 

5 

(b) 

(2) 

6 

(c) 

(2) 

7 

(c) 

(3) 
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£ r ft/s 

2.5 3.5 ' * ’ 7.5 

z R m/s 

Figure 21. Comparison of ^(Zr) for Alternative Lull Pattern 
Directed CROD Letdown Guidance Scheme. 
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|y| 

ft/s 

) x 100% 

Letdown Guidance Strategy it 

0 

1 

2 

3 

B 

5 

6 | 7 

i 

2 

54 

52 

50 

52 

52 

48 

51 

43 

4 

19 

18 

16 

17 

17 

•4 

16 

11 

5 

9 

9 

■ 

9 

9 

7 

8 

4 

7 

2 

1 

1 

1 


1 

1 

0 


TABLE 4 
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TABLE 5 



JF(Z R > x 100% 

*R 



Letdown Guidance Strategy 

0 


i't/s 

D 

fl 

2 

3 

4 

5 

6 

7 

8 

22 

22 

21 

22 

19 

12 

11 

10 

10 

9 

9 

8 

9 

6 

2.5 

1 

0 

12 

3 

3 

3 

3 

1 

0 

H 

■ 


49 





























the CROD letdown profile le resumed to touchdown. In option (3), any 
ongoing CROD profiles are aborted at the instant of lull termination. 
Survival functions t^z^) tor letdown guidance strategies 04, #5, 

#6, and #7 sre shown in Figure 21. These plots show that touchdown 
vertical velocities are lowest for option #7 and monotonically increase 
for option #6, #5, #4 and 0, in that order. 

• ^ 

Survival functions for | yj and j$[ at touchdown are given in 
Tables 3 and 4 respectively. Table 5 presents the results of Figures 
20-21 for^T(zg) in a tabular format, results for | y | ) are not in- 
cluded because there were not any noticeable differences in touchdown 
performance across the? eight letdown guidance strategies. A review of 
these data indicates that the most significant improvements in touchrown 
performance occur for letdown guidance strategy #7. For example, s.rategy 
07 decreases the probability of exceeding touchdown vertical impact 
velocities of 10 ft/s from 9 % for the CROD letdown strategy (#0) to '.OX, 
an order of magnitude improvement. Similar improvements are obtained in 
touchdown values for relative lateral impact velocity |yj and relative 
roll rate | J>| . 

In summary, the results clearly demonstrate the inherent potential 
of a go/no-go letdown guidance strategy which only allows letdown pro- 
files that begin and terminate during lull periods. Fine tuning of such 
a guidance methodology would be necessary with the actual shipboard 
landing scenario under consideration. 
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5. CONCLUSIONS AND RECOMMENDATIONS 


A computer simulation of the VTOL shipboard landing task was 

developed and used to investigate alternative letdown guidance system 

concepts. The NAVAIRENGCEN Ship Motion Computer Program [1] was used 

to generate simulated ship/deck motion time history data for a Spruance 

class destroyer (DD963) operating in rough seas (sea state > 5), at 

a speed (V ) of 25 kts with a heading Cj ) of 120°. The program uses 
s s 

a sum of sinusoids approximation to the ship motion spectra in each of 
the six degrees of freedom for generating the time response data. 

Touchdown performance was evaluated in terms of the survival 
functions (i*e.. Probability of variable X exceeding a level x or 
Pr [X > xj) for four key variables - vertical impact velocity z R (s z p - z ) , 
magnitudes of lateral impact velocity (|y|) and acceleration (Jyj), and 
magnitude of relative roll rate |$[. The touchdown performance of a 
constant rate-of-descent (CROD) letdown scheme was used as a baseline 
against which the performance of alternative letdown guidance systems 
was compared. 

Two categories of alternative letdown guidance concepts were con- 
sidered: (1) designs that exploit the existence of quiescent ship motion 
called lulls, and (2) designs that do not depend upon a lull/swell 
pattern in ship motion. The latter approach requires the development of 
letdown control laws based upon ship motion following and/or forecasting 
schemes. However, actual ship motion time histories exhibit a clearly 
recognizable pattern of alternating lulls and swells, which may be used 
for improving touchdown performance. Therefore, this study was restricted 
to the investigation of alternative ship motion pattern directed 
letdown guidance algorithms. A real time algorithm for lull/swell 
classification based upon ship motion pattern features (e.g., phase- 
plane method) was developed. The classification algorithm was used to 
command a go/no go signal to direct the initiation/termination of an 
acceptable landing window. 
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The temporal characteristics of the simulated ship motion time 
histories were shown to be intimately dependent upon the number of 
sinusoids (N) used to represent the motion spectra. In particular, 
the lull/swell pattern of ship motion becomes increasingly evident with 
the growth in the number of sinusoidal components. The results show 
that a sum of six sinusoids {as suggested in reference 5] representation 
does not display a pattern of alternating lull/swell intervals. There- 
fore, a sum of seventy sinusoids was selected as being adequate for 
describing the known ship motion characteristics. 

Results of the shipboard landing simulation show that substantial 
improvements over CROC touchdown performance may be achieved by using 
a lull/swell pattern directed letdown guidance where touchdowns are 
permitted only during lull periods. For example, the probability of 
exceeding touchdown vertical impact velocities of 10 ft/s is reduced 
from 9% for the CROD strategy to 0% for the lull/swell pattern directed 
guidance strategy - an order of magnitude improvement. Similar improve- 
ments are obtained across a wide range of values for variables z_, |y| 
and | (p j . 


Based upon the results of this study, the following plan for 
further research is recommended: 

1. Analysis of Actual Ship Motion Data. 

The purpose of this task would be to assess the validity of the ship 
motion pattern analysis results obtained from the sum of sinusoids model 
with actual data. Any discrepancies between simulated and actual data 
must be analyzed to determine their impact upon the lull/swell pattern 
statistics and the touchdown performance obtained with the pattern based 
letdown guidance strategy. 

2. Piloted Simulation Experiments . 

The objective of these experiments would be to evaluate the potential 
usefulness of alternative ship motion pattern based letdown guidance 
schemes on manual shipboard landing performance. Four letdown guidance 
options are proposed: 
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a. Letdown with no guidance 

b. Letdown with a go /no go guidance (e.g., ship mounted or 
heads up green/red signal) 

c. Letdown with roll/heave stabilized horizon bar display 
(ship mounted or heads up display) 

d. Letdown with a go/no go guidance and roll/heave stabilized 
horizon display (i.e., options b and c) . 

The purpose of the option (c) would be to determine if the pilot is 
capable of learning the lull/swell pattern from observing the ship 
motion with respect to the roll/heave stabilized horizon bar. 

3. Analysis of Data from Piloted Simulation Experiments . 

The objective of this task would be to analyze the data gathered 
from the piloted simulation experiments with a view towards determining 
the relative merits of the four alternative letdown guidance laws. An 
understanding of this data would be used to develop an improved letdown 
guidance scheme. 
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